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Abstract 

We calculate the quantum tunneling rate of a ferromagnetic particle of 
~ lOOA diameter in a magnetic field of arbitrary angle. We consider the mag- 
netocrystalline anisotropy with the biaxial symmetry and that with the tetrag- 
onal symmetry. Using the spin-coherent-state path integral, we obtain approx- 
imate analytic formulas of the tunneling rates in the small e(= l — H/Hc)-liniit 
for the magnetic field normal to the easy axis {6h = it/2), for the field oppo- 
site to the initial easy axis {6h = tt), and for the field at an angle between 
these two orientations {it/2 <^ Oh tt)- In addition, we obtain numerically 
the tunneling rates for the biaxial symmetry in the full range of the angle 6h 
of the magnetic field {it/2 < 0// < vr), for the values of e = 0.01 and 0.001. 
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I. INTRODUCTION 



A quest to understand the macorscopic quantum tunneling of ferromagnetic particles has 
been an important topical issue of intensive theoretical and experimental studies. 
Since the magnitude of the total magnetization M is frozen up at sufficiently low temper- 
ature, the direction M of the total magnetization becomes the only dynamical variable. In 
the absence of an external magnetic field, this direction is subject to the magnetocrystalline 
anisotropy which possesses its easy directions depending on the crystal symmetry. In this 
situation, M has at least two equivalent stable orientations which leads to the oscillation of 
M between them. This phenomenon is called macroscopic quantum coherence(MQC). |l|]- 
1^ Simple analysis shows that in this case the height of barrier is too high to observe MQC 
in experiments. However, by applying a magnetic field H along the direction which lies 
halfway between two equivalent stable orientations, we can control the height and the width 
of barrier and make MQC observable. On the other hand, if we apply a magnetic field in a 
proper direction, the symmetry of two easy directions is broken. Then, one of the two stable 
orientations becomes metastable. In this situation we can obtain the optimal condition for 
the observation of tunneling, which is called macroscopic quantum tunneling(MQT). 
The ferromagnetic particle is typically a single domain with as many as 10^ — 10® magnetic 
moments, which is a macroscopic number of particles. During the dynamical process we 
can not neglect the infiuence of the environment such as phonons, [Q nuclear spins, |^ and 
Stoner excitations and eddy currents in metallic magnets. Even though some of these 
couplings are noticable, it has been reported @]- that they are not strong enough to make 
MQC or MQT unobservable. A few experimental attempts have been made to observe 
the MQC or MQT of a large single ferromagnetic particle or a collection of magnetic par- 
ticles. At present, it is not easy to perform a direct comparison between theoretical and 
experimental results because of stochastic behaviors of the systems. 

Recently, Zaslavskii studied the uniaxial anisotropy which is the simplest case, mapped 
the spin problem onto a one- dimensional particle system, and then obtained the tunneling 
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exponent, the preexponential factor and their temperature dependences in the low barrier 
hmit with a magnetic field applied at some angle to the anisotropy axis. Later, Miguel 
and Chudnovksy |Q| performed the calculation based on the imaginary time path integral 
method, and demonstrated that their result of the tunneling exponent in MQT coincides 
with the Zaslavskii's result in the uniaxial symmetry, for its dependence on the direction 
and magnitude of the applied magnetic field. They discussed the tunneling rates at finite 
temperature and suggested the experimental procedures whose results can be compared with 
their theorectical results. In this paper, we extend the above calculations to the biaxial and 
tetragonal symmetries by applying the instanton technique. Since the result of the biaxial 
symmetry is a generalization of that of the uniaxial symmetry studied by Zaslavskii, and 
by Miguel and Chudnovsky, we can compare our results with theirs by taking the vanishing 
limit of the transverse anisotropy constant. Also, we will explain that MQC and MQT can 
be consecutively observed by changing the direction of if, and discuss their dependence on 
the direction and the magnitude of H. 

This paper is organized as follows. In Sec. |T|, we briefiy discuss the theory of MQT and 
MQC in a ferromagnetic particle based on the standard instanton method. In Sec. ^ and 
|V], we consider the tunneling rate for biaxial and tetragonal symmetry in a magnetic field 
with a range of angles ti/2 < 9h < vr. We derive approximate formulas of the tunneling 
rates in the small e-limit for three angle ranges {9h = 7r/2, -k /2 <^ 9^ <^ -k, 9h = t^)-, and 
present the ^i^-dependence of the WKB exponent in the full range of angles (7r/2 < 9^ < tt). 
The conclusions are given in Sec. 0. 

II. THE MQT AND MQC IN A FERROMAGNETIC SYSTEM 

The tunneling rate F of a ferromagnetic particle escaping from a metastable state in 
MQT or the oscillation rate A between double wells in MQC has the relation 

(r or A) oc exp(-fi). (1) 
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The WKB exponent B in Eq. (|l]) is approximately given by U/hub, where U is the height of 
barrier. The barrier frequency is the frequency of small oscillations around the minimum 
of the inverted potential, and characterizes the width of the barrier hindering the decay 
process, uj^ is inversely proportional to the effective mass of the magnetic particle, where 
the mass is induced by the transverse component of the magnetic field or by the transverse 
magnetic anisotropy constant. Therefore, it is necessary that the magnetic particle has a 
transverse magnetic anisotropy or we apply a transverse magnetic field, in order that MQT 
or MQC of the magnetic particle is possible. For calculating F or A in Eq. ([l|), we introduce 
the angles 6 and for the direction of M in the spherical coordinate system and employ the 
instanton method, in which the Euclidean action is written as 



S[9{r), 0(r)] = V cos d)'^ + E{d{T), 0(r))]rfr, (2) 



where V is the volume of the magnetic particle, 7 = gfis/h, fiB the Bohr magneton and 
Mq the magnitude of the magnetization. The first term in Eq. is the topological Wess- 



Zumino term and the second term is the energy density which is composed of the 
magnetic anisotropy energy and the potential energy given by external magnetic field. 

The action produces from 6S = the classical equation of motion for M called the 
Gilbert equation [|rT|, which in the Euclidean space is written as 

.dM dE 

I— = --fM X — . (3) 
dT ' dM 

We note that the action (0) describes the 1 © 1 dimensional dynamics in the Hamiltonian 



formulation, which consists of the canonical coordinates and = 1 — cos^. [|12| , p!3| 

In the standard instanton method, the tunneling or the oscillation rate is given by the 
formula 
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^xp(-|^), (4) 



where Sd is the classical action evaluated by using the solution of the equation of motion 
derived from 5S = 0. In Eq. (W) Cq is the preexponential factor which is originated from 



the S'^S term, ut a precession frequency in MQT and ujc an oscillation frequency in MQC, 
where cut and lJc are of the order of the barrier frequency Uh- 



III. TUNNELING RATE FOR BIAXIAL SYMMETRY 

In this section we study a system which has the biaxial symmetry, with ±z axes being 
the easy axes of the Hamiltonian. As in Ref. f^, we apply the magnetic field in the xz-plane 
and take the z axis as the initial easy axis when there is no magnetic field. Then the total 
energy E{6, 0) of the system is given by 

E{e, (p) = Ki sin^ e + K2 sin^ sin^ 6 - H^Mq sin ^ cos - H^Mq cos 9 + Eq, (5) 

where Ki and K2 are the parallel and transverse anisotropy constants, respectively, and £^0 
is a constant which makes E{9, (p) zero at the initial orientation. As will be shown later, 
the effective mass of the magnetic particle is inversely proportional to a linear combination 
of K2 and transverse magnetic field while there is no exact analog of the kinetic energy 
in the action (|^). Denoting 6h to be the angle between the magnetic field and z-axis, we 
obtain MQC for Oh = and MQT for t{ j1 < Qu < -n with large possible tunneling rates 
depending on the magnitude of the applied magnetic field. 
It is convenient to introduce the dimensionless parameters 

K2 = K2l2K^, H,, = HjHo, H, = H,/Ho, (6) 

where Hq = 2Ki/Mq. Then the total energy (H) is written as 

E{e, 0) = ^ sin^ O + K2 sin^ sin^ - sin 6^ cos - cos 6 + Eq, (7) 

where E{9, 0) = E{9, (f))/2Ki. The plane given by = is the easy plane, on which E{6, 0) 
is given by 

E{e,0) = ^sm^e - Hcos{e -Oh) + Eo, (8) 
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where = HsinOn and = H cos 6h since H is in x2;-plane. We define to be the 
initial angle and 6^ the angle at which the barrier vanishes by the applied critical magnetic 
field. From Eq. we get the following conditions which 6q and 6c should satisfy, 

sin^ocos0o + ^sin(^o-^//) = 0, (9) 
smecCosec + Hcsm{ec-0H) = 0, (10) 
cos(2^c) + He cosiOc - 9h) = 0, (11) 

where He is the dimensionless critical field with which the barrier just vanishes. From Eqs. 
(p!0|) and (p!ID the critical field H^ is given by 

" (sin2/3 0^ + |cos0Hp/=^)3/2' 
From Eqs. (p!0|)- ([T^) 9c is given by 

• 2|cotgg|V3 

""(2'^) = l + |cotg^p/3 - (13) 

-(^^-)-^- l + |cot^.p/3 - (14) 

For the special cases 9h = 7r/2, 37r/4, n, we have 9c = 7r/2, 7r/4, 0, respectively. 

The practically interesting situation is when the barrier height is low and the width is 
narrow in order to have the tunneling rate large. Such a situation is realized when the value 
of e = 1 — i//ifc is small. For the small value of e, Eq. @ becomes 

sin(2^e) ie-^r]^)-r] cos(2^,) (2e - r]^) = 0, (15) 

where i] = 9c — 9q which is expected to be small for e <C 1. As will be seen later, detailed 
calculations show that t] is of the order of ^/e. Then the order of magnitude of the second 
term in Eq. (|T5|) is smaller than that of the first term by y/e and the value of rj is determined 



by the first term to have rj ~ y 2e/3, except for the 9^ values near 9^ = tt/2 and vr. However, 
when 9h is very close to 7r/2 or tt, sin(2^c) becomes almost zero as shown in Fig. 1, and the 
first term is much smaller than the second term in Eq. (pTSp. Then the value of t] is obtained 
from the second term when ^ 7r/2 or vr, and r] is given by r/ ^ -\/2e for 9^ ~ tt/2 and 
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^ for 9h — from the detailed analysis of Eqs. (|9|)-([TT|). After a little manipulation, 
we obtain from Eq. (|^) an approximate formula of E[d, 0) in the limit of small e given by 

E{5, 0) = i^ssin^ 0sin2(^o + 5) + i?xSin(^o + - cos0) + Ei{5), (16) 

where E{6, 0) is written as E{5, 0) by introducing a small variable 5 = — and Ei{5) is 
a function of only 5 given by 

E,{5) = \ sin(20e)(35^r^ " + ^ cos(2^^,) [^^(e - ^r/^) + - j]- (17) 



As previously discussed for rj, even though the cos(2^^c) term in Eq. {\[7\} looks smaller by a 
factor of 7] which is of the order of ^/e, the second term can not be neglected near 6h = 7r/2 
and TT because sin(2^c) is almost zero for these regions of 6h- When we assume that |0| <^ 1, 
from the energy conservation E{6, 0) = in Eq. (0), is approximately given by 

Since 5 is of the order of -\/e as can be seen from Eq. (p^Tj), Ei{5) is 0(e^/^) or less. Therefore, 
is very small from Eq. (|18|) whose validity is for the full range of angles 7r/2 < ^^j^ < tt in 
biaxial symmetry. However, in case of tetragonal symmetry the magnitude of is not small 
about 6*// = vr, as can be seen in Eq. (|75D. In such a situation it is not possible to expand 
E{5, 0) as powers of and to reduce the classical equation to the one-dimensional equation 
like Eq. (||). 

Since we confirmed that |0| is very small in case of the biaxial symmetry which we study 
in this section, we have the approximate formula of E{5, 0) in Eq. ([16|) given by 

E{5, 0) = [K2 sin2(^o + 5) + sin(^o + 5)]0' + E,{5). (19) 

By introducing a new scaled time variable f = ujqt with ujq = 2'jKi/Mo, the Euclidean 
action (^ becomes 

S[6{f), 0(f)] = hjj - cos(^o + 5)^] + E{6{f), 0(r))}rff , (20) 



where J = VMo/h'y. In the following we also use the scaled angular frequencies Ui defined 
by uJi = Ui/uJo {i = b, c, t). 

There are three approaches for the calculation of Sd from the action (pOl). Firstly, if |0| 
is small, we can perform a Gaussian integration over the variable in the path integral and 
reduce the system to that with only one variable 6. Then it is possible to do the rest of the 
calculation by using the standard instanton method. However, if |0| is not small like in 
the tetragonal symmetry case with Oh = t^, the integrand of Eq. (^) can not be reduced to 
the form which allows the Gaussian integration over 0. Secondly, we can use the fact that 
E{6ci, (pel) = and that {d(f)/d6)d6 can be substituted for {d(f)/df)df in the Euclidean action 
(pop. Here, we do not need to know the explicit solutions of the classical paths Sci{f) and 
0c«(^) ^he classical action Sd because we can obtain d(f)/d6 = {dcj) / df) / {d5 / df) from the 
Euler-Lagrange equations given by Eq. (|21|) and ( ]22|) below. However, this approach can 
not give the value of the preexponential factor, for which we need the explicit solution of the 
bounce in MQC or the instanton in MQT as a function of f. [|14| Thirdly, we directly solve 
the coupled equations of motion for 5 and given by Eqs. ( PT| ) and (^) by incorporating 
E{Sd, (pd) = 0. Even though it is hard and sometimes rather tedious to obtain the solution 
as a function of r from Eqs. ( PT| ) and (^), and E{6d,(pd) = 0, this methods provides a 
complete result for the tunneling rate F in MQT or the oscillation frequency A in MQC. In 
this approach, the preexponential factor can be obtained from the second variational term 
of the action in the spin-coherent-path integral, and the classical action Sd is obtained 
by the direct integration of Eq. ( pOD over r with the explicit solution 6d{f) and 0c«(^)- While 
the third method is the most general approach to obtain the tunneling rate, it is reduced 
to the first one if |0| is small enough. Also, the second method is useful for checking the 
result of the first or the third method for the WKB exponent. In biaxial symmetry, due to 
101 <^ 1 we apply the first method to obtain the explicit instanton or bounce solution as a 
function of f, which is used to calculate the WKB exponent and the preexponential factor. 
On the other hand, since the magnitude of is not small in tetragonal symmetry case with 
6h = vr, we need to make use of the third method for the calculation of F. In addition, by 
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using the second method we check the results of the WKB exponent from the first or third 
method for each case. 

The classical trajectory of 6 and (p is determined by the Euler-Lagrange equation derived 
from the action pO|), 

From Eqs. (|19D, ( ^T]) and (|2^) , the equation of motion for the instanton 6ci{f) in the biaxial 
symmetry case becomes 



where the effective mass is given by 

• ^ +2i^2]'^. (24) 

sm[UQ + del) 

Now let us study the tunneling rates at three different magnetic field directions of = 71/2, 
TT, and tt/2 <^9h <^n. 



A. 9h = it/2 



For 6h = 7r/2, 6c = 7^/2 from Eqs. ( [TB| ) and ([T^), and rj = v2e as explained below Eq. 
(|T5|). Then, from Eqs. (0) and (|2^) we get the approximate forms of the effective potential 
energy Ei{6) and the effective mass M in the reduced one dimension as 

= l^'i^^ - v^)^ (25) 
M ={H, + 2K2)-\ (26) 

where Hx ~ He- Here we notice that MQC can not be generated without an external field 
along the x-axis or a transverse anisotropy constant i^2- Eq- (^) with Eqs. and 
has the solution called the bounce which is given by 



Sci{f) = V2e[l + tanh(cjer)], (27) 



where Uc = y e/2yl + K2/ Ki. Then from Eq. ([T8|) we get the classical path for (p given by 

^ci (r) = -^e^=i== [ J, A • (28) 
^1 + ^ cosh [uJct) 

We can calculate the classical action by integrating Eq. (^) with the above classical tra- 
jectory, and the result is given by 

4a/2 6^/2 

S^i = ihJ)— ; ■ (29) 

From Eq. (p5|) and Fig. 2, we get the height of barrier as Ei{= U/2KiV) = e^/2 at 
5m = V^e and the oscillation frequency around the minimum of the inverted potential 



—Ei{6) as ujb{= y —Ei{6m)/M) = y/2ujc- Then, as mentioned in Sec. 0, the WKB exponent 
B{= Sci/h) is approximately given by 

Br^^ = - ' (30) 



which agrees up to the numerical factor with the result in Eq. (pOD obtained by using the 
explicit instanton solution. 



B. it/2 «: 61^ «; vr 

For 7r/2 ^ ^ 7r, the critical angle 9c is in the range of ^ ^ 7r/2 and t] ~ ^2e/3. 
Then from Eqs. (|T3|) and (|T^) we get 

Ei{6) = ^sm{2ec){V6'e6' - 5=^), (31) 

M =[^ + 2K2]-\ (32) 
sm 6 c 

where = He sin 9h- The classical equation of motion Eq. (pSf) gives the instanton solution 

as 
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cosh^ (uJtT) 



where Eq. ( p!8D has been used to get 0c/ (t), and the dimensionless frequency Ut is given by 





cot 6h 


|l/6 


1 + 


\cote 


'h\ 


2/3 



Then the classical action is found to be 



= (36) 



It is noted that Eq. (|36| ) without agrees with the classical action in Ref. |Q| which 
studied the uniaxial symmetry case. By using Ei{5m = 2v^/3) = sin(20J(6e)3/V27 and 



Ub = y —E"{6m)/M = 2u!t given from Eq. (^Tj) and Fig. 2, we approximately obtain B as 



which is consistent with Eq. (|3^) up to the numerical factor. 



C. eH = TT 



In case of 9h = vr, we had 9c = and = 0. Then we have 



EiiS) = lies' ~^-^), (38) 
M = {2K2)'' (39) 



which gives the classical path obtained from Eq. (^) as 



Ur) = If , , (40) 
cosh(c<;tr) 

4>d{r) = -iJ tanh(u;tf ) + rni, (41) 
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where n = or 1, and cut = 
conservation of energy E{6ci, (pd 



K2e/Ki. Here we note that (pdiT) is obtained from the 
: 0. Then the corresponding classical action is given by 



Sci = (hJ) 




8 iKi 



We note again that 



B 



U 



(42) 



(43) 



htOb 4 

which is obtained from the fact that E{5m = "\/2e) = e^/2 and ujb = V^^t, is consistent with 



Sci in Eq. (42) up to the numerical factor. 



D. Prefactor and Discussion 

In order to complete our study of MQT and MQC in the biaxial symmetry, we need to 
calculate the preexponential factor in Eq. (^). Since in the biaxial symmetry the problem 
can be reduced to the one-dimensional one due to the smallness of \(f)\, the calculation of the 
preexponential factor can be performed by using the well-known instanton method developed 
by many authors. |T^- |]TB[ However, in many other symmetries it is not possible to reduce 



the problem to the one-dimensional one by directly integrating over one of two variables 
such as (j), because the magnitude of (j) is not small enough for the action to be expaneded as 
powers of (p for Gaussian integration. This situation will be seen in the tetragonal symmetry 
case with 6h = tt. In order to treat such a case, Garg and Kim Ijl^ studied the formalism for 
evaluation of the prefactor based on the spin-coherent-state path integral. Here we explain 
briefly the basic idea of this study. Expanding the action S[9{t), (p^r)] in terms of 9i and 
01 about the classical path, where 6{t) = 9^ + 9i and 0(r) = (p^L + 01? we obtain the action 
as S[9, (j)] ^ Sci + ^'^S with being a functional of 9i and 0i. Then Gaussian integration 
can be performed over 0i, and the 9i path integral is now cast into the standard form of 
one-dimensional potential problem. What we need in this calculation is the approximate 
form of d9ci/dT for large r. 
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In the biaxial symmetry, by using Eq. (^) for 6h = vr/2, Eq. (|33|) for 7r/2 ^ 6*// ^ vr, 



and Eq. (EDI) for 6^ = tc, we obtain the complete anlaytic forms of the tunneling rate for 



MQT or the oscillation rate for MQC, and the results are summarized in Table 11(a). p!8|JT9 
We note that in Table 11(a) the preexponential factor is multiplied by a factor of two in 
case of 6h = because there exist two classical paths which give the same action. Also it 
is noted that the e behavior of the WKB exponent B is given by e^^^ for 9h = 7r/2, e^^^ for 
7r/2 <^ 9h ^ vr, and e^^"^ for 9h = tt. It is seen by taking K2 to be zero in Table II ||2^ 
that the tunneling rate in the uniaxial symmetry vanishes for Oh = vr. This situation can 
be understood from the fact that for Oh = the Possion bracket {p^, if} which determines 
the dynamics of the spin system with the Lagrangian (H) is zero because in this case the 
Hamiltonian H becomes a function of only p^, which gives dp^/dr = d{l — cos 9) /dr = 



and then 6 is constant in time. 13 



We obtained the instanton solutions for the full range of 9h (vr/2 < 6^^ < vr) by solving 



numerically the equations of motion ( [2l| ) and ( ^2]) or equivalently Eq. (p3D. Then by using 
the obtained instanton solutions we calculated the classical action from Eq. ( PPD to obtain 
the WKB exponent B{= Sd/h). In Fig. 3 we present the instanton solutions with e = 0.001 
and Ki = K2 for several values of Oh which we obtained by numerical calculations. We 
also obtained the ^^-dependence of B with e = 0.01 and e = 0.001 for 7r/2 < ^^j^ < tt, and 
present the result in Fig. 4. As is noted in the figure, the maximal value of B is at about 
6h = 2.78(~ 159°) and the approximate analytic result obtained in Eq. ( p6D is almost valid 
in the range of angles 120° <6h< 170°. 



IV. TUNNELING RATE FOR TETRAGONAL SYMMETRY 

In this section we study the tetragonal symmetry whose anisotropy energy is given by 
Ea{e, 0) = Ki sin^ e + K2 sin^ 6 ~ K'^^ cos(40) sin^ 9, (44) 
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where we once again take the easy axis to be ±z for Ki > 0. When we apply H in the 
xz-plane as in the previous section, the total energy becomes 

E{9, 0) = Eaie, 0) - H^Mo sin ^ cos - H^Mo cos 9 + Eq, (45) 

where we assume that Ki > and \K2 — K',^ <^ 1. ||2l| Here we also use the dimensionless 
parameters defined in Eq. (^). By choosing K'2 > 0, we take = to be an easy plane in 
our calculations. In the = plane the scaled total energy is written as 

E{9, 0) = - sin^ e + {K2- R'^) sin^ d-H cos{9 - 9h) + Eq, (46) 

where = K^/^Ki. The initial angle 9q and the critical angle 9c defined in the previous 
section are determined by the equations 

sin 9q cos 9q + H sin(6'o — 9h) + 4(i?2 - ^2) sin' 9o cos 9o = 0, (47) 
sin 9c cos 9c + He sin(^, - 9h) + 4(7^2 - K) sin' 9c cos 9c = 0, (48) 
cos{29c) + He cos{9c - ^h) + 4(ir2 - iT^) (3 sin^ 9c cos^ 9c - sin^ 9c) = 0. (49) 



Using Eqs. (^) and (^5]), we obtain the critical magnetic field as ||2^ 

1 . 4(i?2 - K2 



Hr 



\ . 2/3 ]- (50) 
1 + cost^H-r' 



(sin2/3^^ + Icos^nP/')'/' 1 + Icos^i^p/sJ 
In the small e-limit (e = 1 — H/Hc), by using Eqs. (^) and (^9|) we obtain the approximate 
equation for ?](= 9c — 9q) given by 

- e[2^,sin(^^, - 9h)] + r^^pi?, sin (0, - ^^^) + 6(7^-2 - sin(4^,)] 

+r/{e[2i?,cos(^e - 9h)] - v^[HcCos{9c - 9h) + 8(7^2 - cos(4^e)]} = 0. (51) 



The order of magnitude of the first two terms in Eq. (|5TD is higher than that of the last 
third term by 0{y/e). Therefore, for 7r/2 ^ ^^j/ ^ vr we obtain rj from the first two terms of 
Eq. (|1]) as 

r/ = y^[l + 4(K2-i^;)cos(2^,)], (52) 
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and dc from Eqs. (||) and (H) as 



. ^ , (53) 

However, for 9^ = 7t/2 and vr, since sin(6'c — Oh) and sin(4^c) are equal to zero, the first two 
terms in Eq. (^l]) vanish. In these cases the value of rj is determined by the last third term, 
which gives 

0c = V^V2'e[l-A{K2-K',)] for = |, (54) 
= 0, ri = for Oh = vr. (55) 

For reference, we note that t] is approximately given by 

for the vaule of 9h around 9h = 7r/2. 

In the small e-limit the approximate form of E{6, 0) in Eq. (|46|) is written as 

- ^[H, cos{e, - Oh) + 8{K2 - R',) cos(40,)] [6'{e - ^v') + V^' - j] 

+ 4(^2 - R2) cos{Ae,)6^e, (57) 

where 6 = 6 — 6q which is small in the small e-limit. Then the total energy becomes 

E{6, 0) = R'^ll - cos(40)] sin^(^o + + H^{1 - cos0) sin(^o + 5) + Ei{6), (58) 

which will be discussed for three angle ranges, 6^ = n /2, n and ti /2 <^ 6^ <^ n. 



A. Oh = vr/2 



Using the equation of motion for the classical trajectories, Eqs. 
classical path at 6*// = 7r/2 is given by 



and (|22|), the 
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5,i{f) = V2t[l - 2^^^— ^][1 + tanh(cD,f)], (59) 

Mr) = -^e[l - 2 ^^^^^^^ ][ , _J , (60) 

Ki cosh [uJcT) 

where Uc = \J e/2[\ + 2[K2 + K'2)/Ki\. In this case the approximate form of Ei[5) and the 
effective mass become 

[i + i2(i?2-i^;)],2 



Ei{S) = Y - 2v/2e[l - A{K, - K,)]}\ (61) 

M = ^ ^ — —. (62) 

ifc + I6K2 ^ ' 

Since 1 0^(^)1 ^ 1, we can obtain a reduced one- dimensional action by performing the 
Gaussian integration over 0. In this case Eq. (^) is the effective potential energy in the 
equation ( p3l) for the bounce 5ci{f). Using Eqs. (|^) and (|60|), the corresponding classical 
action is given by 

S^, = {hJ)\V2e"\l~A^). (63) 

From Eq. (|6lD and Fig. 2, the height of barrier and the barrier frequency for 9^ = 7r/2 are 
given by 

Ei{5.m) = 2e'^[l - A{K2 - K'^)] and cUb = V2lJc, (64) 

where 5m = v^[l ~ 4(^2 ^ -^2)]- Then the approximate form of the WKB exponent 
becomes 



which up to the numerical factor is consistent with the action Eq. ( pB| ) obtained by the 
calculation with the explicit bounce solution. 



B. Oh Oh <^ TT 



In this case the potential energy and the effective mass in the reduced one dimension are 
approximately given by 
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The classical trajectory is given by 



Sciif)= (68) 
cosh [uotT) 



.(6e)3/1cot^.|V^^^^^ 

cosh [UJtT) 



K'^ 2,2-3|cot0iyp/3 i;^2-i^2 



1 + I ootfa|V3|l - 411^ + )(^)|. (69) 

where the precession frequency cjf in MQT is given by 

The corresponding classical action becomes 

[^~^K: + 3^1 + |cot^Hr^^~^^^- 
Using the barrier height Ei{5m = 2v^/3) = 4a(6e)3/^/27 where a = — [ffcsin(6'c — + 
2{K2 — K'2) sin(46'c)]/2, and the oscillation frequency Cub = 2iUt, we obtain the approximate 
form of the WKB exponent as 

which is consistent with Eq. ( |71| ) up to the numerical factor. 



For 9h = the total energy is written as 

= K',[l - cos(40)]5^ + l{e52 - ^[1 - 8(i^2 - K',)5']}. 



(73) 



From the Euler-Lagrange equation (^Tj) and ( P^ for 5(r) and 0(t), the classical trajectory 
is given by |]17 
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\\k + K'^ cosh(4a;tf ) ' 

n 

-luJtT + -n, 



(74) 
(75) 



where n = 0,1,2,3, K = Ki/A — K2, and ut = e- Eq. (^) shows that in the present case 
101 ^ 1 is not vahd and that we cannot expand E{6,(f)) as powers of (p, which means that 
we cannot reduce the problem to the one- dimensional one like Eq. (|23|) . Thus the classical 
action should be obtained directly from the solution ([7^ ) and (^) of the Euler-Lagrange 
equation, which becomes 



4 Ka K2 



(76) 



where Ka = y — {K'2)^. Even though the effective potential energy and the mass in 
one-dimensional form are not appropriate for the present case, the dependence of i? on e can 
be derived from the values of Ei{6m) and ut, in Table I, which gives e^~^ from ^i((5m)(oc e^) 
and LjJoc e). 



D. Prefactor and Discussion 

Following the calculations used in the previous section, we can obtain the preexponential 
factor Cq in Eq. (^). As can be shown in Table II, for 6h = v/2 and 7r/2 ^ ^ vr the 
values of Co in the tetragonal symmetry coincide with those in the biaxial symmetry. This 
situation is understood by the fact that the prefactor Cq is determined only by the shape 
of the effective potential energy in the reduced one dimension, and the shapes are the same 
for the two symmetries as can be seen in Eqs. (p5D , (|3T|) , (|6TD, and (|66D . However, the 



situation is different for Oh = because |0| <^ 1 is not valid. Thus, the factor Co should 
be determined from by using the approximate forms of 5ci, Es^, E^^, Ess, and so on, in 
the tetragonal symmetry case with 9h = tt, which leads to the result of Co in Table 11(b). 
It is possible to perform the same numerical calculation for given values of the parameters 



K2 and K'2 in the tetragonal symmetry in the same way as we did in the biaxial symmetry 
for the full range 7r/2 < 6h < vr. The results of the detailed numerical calculation for the 
dependence of B on Oh will be presented elsewhere. 

V. CONCLUSIONS 

In summary we obtained the approximate analytic forms of the oscillation rate (MQC) 
for 6h = 7^/2 and of the tunneling rate (MQT) for 7r/2 ^ ^ tt, and = vr in the biaxial 
and tetragonal symmetries. Also, the ^//-dependence of the WKB exponent B for the full 
range 7t/2 < 6h ^ n with the values of e = 0.01 and 0.001 was calculated numerically in 
the biaxial symmetry. For future works it will be needed to compare the theoretical results 
of the ^//-dependence of B for given e and those of the e-dependence of B for given Oh 
obtained in this paper for the biaxial and teteragonal symmetries, with the experimental 
results which can be obtained by the same procedures as those suggested in Ref. p for the 
uniaxial symmetry. 
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TABLES 

TABLE I. The values appearing in Fig 2 for three ranges of Oh (a) in the biaxial symmetry 
and (b) in the tetragonal symmetry. 5m is the angle where the maximum of the barrier locates. 
Here, a = -[HcSm{ec - Oh) + 2^sin(4ec)]/2 where K = {K2 - K'^)/2Ki. 
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TABLE II. The oscillation rate in MQC for 6*// = 7r/2 and the tunneling rate in MQT for 
7r/2 ^ 0// ^ vr and tt (a) in the biaxial symmetry and (b) in the tetragonal symmetry. Cq is 
preexponential factor, ujf, in MQC or ujp in MQT the frequency of small oscillations around the 
minumum of the inverted potential, and B the WKB exponent. The complete form of the rate is 
given in Eq. (§. We define i^;3 = ik+K^)/K'„ K, = {K2-K',)/K,, = (f)'^' (i||gotXp3) 
and Bu = ^^^^ ^ \ cot 0//|^/^e^/^, where Cou and Bu are the dimensionless frequency and the WKB 
exponent for 7r/2 ^ 0// ^ vr in the uniaxial symmetry. [0] 
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FIGURES 

FIG. 1. Comparison of (a) sin(2(9c) with (b) cos(2(9c) in Eqs. and (||). Note that 

sm{29c) = 0, I cos(2^c)| = 1 for both 6h = vr/2 and vr. 

FIG. 2. The plots of the effective potential Ei{5) as a function of 6{= 9 — 6q) for (a) 9h = 7^/2 
(MQC), (b) Sit /A and (c) vr about 6q at which the system is metastable (MQT). The position 6m 
and the height Ei^Sm) of the barrier, and the barrier frequency cD^ are given in Table I for the 
biaxial and tetragoanl symmetries. 

FIG. 3. The ^//-dependence of the instanton solutions in the biaxial symmetry for 
e(= 1 - H/Hc) = 0.001 and Ki = K2. Here, (a) Oh = 90.0002°, (b) 90.01°, (c) 135°, (d) 
179.99°, (e) 180°. 

FIG. 4. The ^//-dependence of the relative WKB exponent B{6h)/B{2>tt/4) in the biaxial sym- 
metry with i^i = K2 for (a) e = 0.01 and (b) e = 0.001. Here, 5(6'// = 37r/4,e = 0.01) = 9.14x10"^ 
and B{eH = 37r/4, e = 0.001) = 5.14 x lO""^. 
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